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A Kerr black hole with mass M and angular momentum J satisfies the extremality inequality J < 
M 2 . In the presence of matter and/or gravitational radiation, this bound needs to be reformulated 
in terms of local measurements of the mass and the angular momentum directly associated with the 
black hole. The isolated and dynamical horizon framework provides such quasi-local characterization 
of black hole mass and angular momentum. With this framework, it is possible in axisymmetry to 
reformulate the extremality limit as |J| < 2M%, with Mh the irreducible mass of the black hole 
computed from its apparent horizon area and J obtained using approximate rotational Killing 
vectors on the apparent horizon. The jj| < 2M% condition is also equivalent to requiring a non- 
negative black hole surface gravity. We present numerical experiments of an accreting black hole 
that temporarily violates this extremality inequality. The initial configuration consists of a single, 
rotating black hole surrounded by a thick, shell cloud of negative energy density. For these numerical 
experiments, we introduce a new matter-without-matter evolution method. 

PACS numbers: 



I. INTRODUCTION 

Kerr spacetimes, representing a spinning black hole 
(BH) in isolation, have a bound on the maximum al- 
lowed angular momentum. If M is the mass of the BH 
and J its angular momentum, the bound or extremality 
condition reads |J| < A/ 2 . An extremal Kerr BH sat- 
urates this condition (i.e. | J| = M 2 ). Kerr spacetimes 
violating the extremality condition (i.e. \J\ > M 2 ) have 
naked singularities instead of BHs. A natural question 
to ask then is whether more general BH spacetimes (e.g. 
BHs in the presence of matter, Maxwell fields, gravita- 
tional radiation and/or other BHs) are also subject to 
extremality conditions, thus precluding the existence of 
super- extremal BHs. 

The existence of extremality conditions in more gen- 
eral BH spacetimes has been investigated in a series of 
papers by Ansorg and collaborators. These studies con- 
sidered axisymmetric and stationary BHs surrounded by 
matter [h, In some cases, Maxwell fields were also in- 
cluded Q ■ Of direct relevance to our work is the study by 
Ansorg and Petroff 0] demonstrating the feasibility of 
configurations with super-extremal BHs. More recently, 
in a study of binary BH initial data with nearly extremal 
spins [f|, initial data with super-extremal BHs were 
constructed; however, the binary black holes (BBHs) 
were enclosed by an apparent horizon (AH) with sub- 
extremal properties, thus effectively acting as a single, 
sub-extremal BH. 

As pointed out by Booth and Fairhurst [f|, at the 
core of any study of BH extremality conditions is char- 
acterizing what one means by mass and angular momen- 
tum associated with the BH, namely distinguishing the 
"local" properties of the BH from those of its environ- 
ment. As shown also in the natural tools for such 
a task are those provided by the isolated and dynami- 



cal horizon framework |7(. Using this framework, Booth 
and Fairhurst Q reformulated extremality conditions in 
terms of restrictions on isolated and dynamical horizons, 
and introduced a parameter that determines how close a 
horizon is to extremality. 

The goal of our work is to investigate the formation 
of a super- extremal BH in a dynamical setup using the 
tools of numerical relativity. Specifically, we consider a 
sub-extrcmal BH surrounded with a spherically symmet- 
ric cloud of negative energy density. As the BH accretes 
the cloud, its mass decreases and its angular momentum 
increases. As a consequence, the dimensionless spin pa- 
rameter of the BH increases. This approach to spinning 
up a BH by decreasing its mass is complementary to the 
more astrophysically relevant case in which a BH spins up 
by gaining angular momentum from the material it swal- 
lows Q . For clouds with enough negative energy density, 
we are able to build super-extremal BHs. The super- 
extremal state is, however, not stable. Non-axisymmetric 
instabilities develop, triggering emission of gravitational 
radiation that carries with it copious amounts of angular 
momentum. 

Our study also introduces a new evolution scheme 
called matter-without-matter (MWM). The MWM ap- 
proach consists of evolving the spacetime geometry using 
the BSSN Q evolution equations without their matter 
source terms, namely the vacuum version of the evolu- 
tion equations. The matter fields are constructed at each 
step from the BSSN constraints (Hamiltonian, momen- 
tum and connection constraints). In addition, the stress 
energy tensor is required to have a form for which the 
matter source terms in the BSSN evolution equations 
vanish. The approach is reminiscent of the constraint- 
violating BH initial data evolved in a previous study [Io| . 
Our MWM method differs, however, from work with 
hydro- without-hydro evolutions [ll|, [l2[ . In those stud- 



2 



ies, the matter hydrodynamics was pre-prescribed and 
was used as input to construct the source terms in the 
evolution equations. 

The paper is organized as follows: In Section UH we 
summarize the extremality conditions and discuss the 
condition used in our numerical experiments to identify 
the emergence of super-extremality. Section lUtl describes 
the MWM evolution approach and the conditions that 
the matter content must satisfy. In Sec. |lVl we describe 
the initial data configurations used in our simulations. 
Results showing the formation of a super-extremal BH 
are given in Sec. |V] In Sec. I VI) we address the de- 
gree to which the null energy condition is satisfied by 
our spacetimes. In Sec. IVIII ,we investigate the late-time 
behavior in terms of the constraints. Conclusions are 
presented in Sec. I Villi The numerical simulations and 
results were obtained with the Maya Code as described 
in Refs. [HI [H, [H| • Subscripts a, b,c, .. will denote space- 
time indices while k, ... will denote spatial indices. We 
use units in which G = c = 1 . 



II. EXTREMALITY CONDITIONS 

There are three notions of extremality in stationary, 
asymptotically flat BH spacetimes (see Q for details): i) 
the angular momentum J and mass M of a BH satisfy 
the bound \ = I J\/M 2 < 1; ii) the surface gravity n of 
a BH satisfies the bound n > 0; and iii) the interior of a 
non-extremal BH contains trapped surfaces, while there 
are no trapped surfaces in the interior of an extremal BH. 

The first notion of extremality (x < 1) compares the 
Christodoulou mass [l5[ M to the angular momentum J 
of the BH. The Christodoulou mass is given by 
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Mh is a local measure of the BH mass, called the ir- 
reducible mass. The mass is computed from Mh = 
y/A/16-n, with A the area of the AH of the BH. The 
angular momentum J is also obtained locally using ap- 
proximate rotational Killing vectors on the AH Q. 
For stationary BHs, the surface gravity is given by Q 



(i-x 2 ) 1/2 



2M H [l + (l-x 2 ) 1/2 } 
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Evident from Eq. ([2]) is the equivalence between the first 
two notions of extremality, namely % < 1 and K > 0. 

The goal of our study is to investigate conditions un- 
der which the extremality conditions X < 1 and k > 
are violated. Using x for this purpose has, however, the 
following drawback. From Eq. <jTJ) , it is not difficult to 
show that 
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Thus, by construction x < 1- It has been then sug- 
gested H, !| that a quantity more suitable to investigate 
extremality is £ = \J\/ (2M H ). In terms of £, the spin 
parameter x an d the surface gravity k take the form 
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respectively. Thus, for £ < 1, one recovers the extremal- 
ity conditions x < 1 an d k > 0. The advantage of the 
new spin parameter £ is that values £ > 1 are allowable. 
Furthermore, for £ > 1, one still has x < 1, but the 
surface gravity, on the other hand, becomes negative. 

In Ref. Q , a new notion of extremality was introduced 
in terms of £: A BH is said to be sub-extremal if £ < 1 
(k > 0), extremal if ( = I (n = 0), and super- extremal 
if £ > I (k < 0). This new notion of extremality is de- 
rived assuming that: i) the spacetime is axisymmetric; 
ii) the null energy condition is satisfied; and iii) the cross 
sections of the horizon are embcddable in Euclidean K 3 . 
A generalization of extremality, relaxing the axisymme- 
try and embeddability assumptions but keeping the null 
energy condition, is also possible [f|. 



III. MATTER- WITHOUT-MATTER 
EVOLUTIONS 

Our study introduces also an evolution scheme that 
we call the matter-without-matter method. Under the 
MWM approach, the geometry of the spacetime is 
evolved using the vacuum or source-free BSSN evolution 
equations. There is no need for matter evolution equa- 
tions. As we shall show next, the matter fields are ob- 
tained at every step of the evolution from the BSSN con- 
straints (Hamiltonian, momentum and connection con- 
straints) together with equations- of- state. To demon- 
strate how the MWM works, we will re-derive the BSSN 
evolution equations, taking explicitly into account a mat- 
ter content that is invisible to these equations. 

The BSSN formulation of the Einstein equations con- 
sists of a set of evolution equations for the conformal 
metric 7^ , the conformal factor </>, the trace of the phys- 
ical extrinsic curvature K, the trace-free part of the 
conformal extrinsic curvature Aij, and the connection 
f 4 = ^ k f ) k = -d 3 f3 (see Ref. [H for details). The 
evolution equations for tfi, 7.^ and K are respectively: 



f 



d a (j) 
d a K -- 



— aK 
6 

-2 a An 



-K- 



(0) 
(7) 

(8) 



V denotes covariant differentiation with respect to the 
physical metric 7^ = e^jij , a is the lapse function, and 
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j3 l is the shift vector; we define d a = dt — Cp, with Lp the 
spatial Lie derivative along (3 l . Above, the source terms 
p and S = 7 y Sij arc obtained from 

p = n a n b T ab (9) 
f = -Y a n b T ab (10) 
Sij = ltl]T ab (11) 

with T a b the stress-energy tensor and with n a the unit, 
time-like normal to the constant t hypcrsurfaces. That 
is, the stress-energy tensor has the following form: 

Tab = P n a n b + 2 j (o n b) + S ab ■ (12) 

Before considering the evolution equations for Aij and 

r l , we need to recall the constraints in the BSSN formu- 
lation. As with any 3+1 formulation of the Einstein field 
equations of general relativity, the BSSN formulation in- 
volves the Hamiltonian and momentum constraints. In 
terms of the BSSN variables, these constraints read: 

e - 5 *V 4 V l e* -^-R+lA ij A ij - ^-K 2 = -2np (13) 
8 8 12 

and 

e-^V^e 6 *^') - \v l K = 8irf , (14) 

respectively. In addition, the BSSN formulation intro- 
duces the following new constraints: 

A = A i i = (15) 
S = 7-1 = (16) 

g l = r + djf j = o . (17) 

Our Maya Code actively imposes the trace-free (fl~5j) and 
unit-determinant (fT6|) constraints. On the other hand, in 
numerical evolutions, the Hamiltonian (|13|) . momentum 
(p~4|) and connection (|17|) constraints are not explicitly 
imposed. One of the main virtues of the BSSN formula- 
tion is precisely that aspect. Namely, in the course of a 
BSSN evolution, the constraints (fl3"]l. (fI3|l and $T7J) are 
preserved within tolerable levels. In our MWM scheme, 
the Hamiltonian, momentum and connection constraints 
play a different role. The Hamiltonian and momentum 
constraints are used to construct p and j l , respectively. 
But more importantly, in MWM simulations (see below) , 
Q % evolves away from Q l ~ and provides a key ingredi- 
ent in determining the dynamics of the matter content. 
For this reason, we need next to re-derive the evolution 
equations for A^ and T l without imposing the Q l = 
constraint. 

For the A^ evolution equation, the only place where 
Q l enters is in the Ricci tensor R^, which is normally 
computed from from Rij = Rij + Rf- with 

Rij = -\l lm did m ^ 

+ l Hl d 3) f k +f fc f fc)fc , (18) 



and 

R% = -2(ViVi0 + 7«V*VV) 

+ 4(v.^V J 0-7 y V fc 0VV) , (19) 

respectively. If Q l ^ 0, the Ricci tensor Rij acquires an 
additional term, namely 

Rij -> Rij = -^l lm did m % 

+ 7 ,m (2T \^j)km + r* m r fe yj 

+ 7fe( ! :^)f fc +f fc f fo - )fc 

- lk{idj)Q k -g k T (m . (20) 

Therefore, Rij is given instead as Rij = R^ + Rfj + 

8 7T Mij with R^ and given by (fTH|> and (fT^|) , respec- 
tively, and 

^^-^(^d ]} G k +g k f {lj)k ) . (21) 

With viewed as an additional source term, one can 
follow the standard BSSN derivation [l|| of the evolution 
equation for A^ and arrive at: 

doAij = e~ A * [-ViVj-a + aiEy - 87ra(5y - M l] )] TF 
+ a(KAij-2A ik A k j), (22) 

where TF denotes the trace-free part of the tensor. 

The remaining equation to consider is the evolution 
equation for T l . Once again, we will derive this equation 
without the assumption that the constraint Q % = holds. 
The starting point is the time derivative of T' = —djj IJ + 

d t r = -djd t f j + d t g l . (23) 

On the other hand, from Eq. ([7]): 

da 13 = 2aA ij +p k d k f j 

- l lk d k p - ^ k d k (3< + ^f 3 d k (3 k , (24) 

which after differentiation reads 
d t djf j = d j {2aA ii ) 

+ /3 k d k djf 3 - dj^dkP + ^dj^d k p k 

- y^djdkP - p ij djdkl3 k . (25) 
Thus, Eq. ((Ml) becomes 

d t r = -2A ij d j a-2ad j A ij 

+ j jk djd k ^ + p ij djd k p k + dtQ 1 . (26) 
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Next wc use the momentum constraint ([14]) to eliminate 
the term involving djA^ , and rewrite (]2T5)) as: 



-2A ij d 3 a - l&iraf j j 3 



-2a \-r l jk A^ - QA^drf + ^f 3 d 3 K 

+ ^ k d 3 d k p l + \f ] d 3 d k p k + d t g l . (27) 



the Hamiltonian (jT3|) . momentum ([14)) and connection 
(fTT]) constraints. Since the MWM method can be also 
viewed as evolving constraint violating data, there are 
no guarantees that the method is capable of producing 
stable evolutions. A general proof of the conditions that 
the initial data have to satisfy in order for the MWM to 
yield stable evolutions is beyond the scope of this study. 
We have found, however, that the set of configurations 
and evolutions for the present study on super-extremality 
were all numerically stable. 



With the help of d^ ij = -P + Q\ we obtain 
d t f* = 2. 1"/;,,. Hirrn- '■'./, 



2 a -T) k A> k - 6A v d j $+- s y ,3 d j K 



+ pdjT* - Y'<), f + - I'Oj i 1 

+ ^djdk? + p ij djd k p k + d t g l , (28) 



which can be rewritten as 



IV. INITIAL SETUP 

Our study consists of a single, spinning BH punc- 
ture [l7T | enclosed by a thick, spherically symmetric, 
shell cloud of (in most cases negative) energy density. 
The shell surrounding the BH has a Gaussian profile 
and is initially static, with stress-energy tensor given by 
Tab = p(n a nb — 7ab/3). Specifically, with this choice of 
stress-energy tensor, initially f = Q l = d Q l = 0. 

Vanishing j l implies that the momentum constraint 
reduces to the vacuum case. One can thus directly use 
the Bowen and York extrinsic curvature solutions for a 
spinning puncture (l8j : 



d n r = -2A v d, 



-2a ( -V) k A ik - (i.V'H,o + --"«;, A 



(29) 



with both r* and Q l treated as vector densities of weight 
2/3. Equations @, ©, ©, fFfy, 0, (EH and 0® 
constitute the basis of our MWM method. 

From the form of the stress-energy tensor (|T2"j) . we see 
that without any further assumptions, our set of equa- 
tions is not able to determine SV? . The essence of a MWM 
evolution is then to impose that matter evolves in such 
a way that the source terms in © , (|2"2"j) and (|2"9"j) vanish, 
namely 



P + S = 



TF 



(Sij - My) 

d o g l - 16 7rae~ 4 V = 0. 



(30) 
(31) 
(32) 



These conditions imply that the stress-energy tensor (TP2")) 
takes the following form: 

Tab = P (n a n b - ^lab^j + 2j (a n b ) + 7o7&My F • (33) 

The obvious advantage of the MWM method is the 
direct use of a vacuum (e.g. black hole) BSSN evolution 
code to construct the spacetimc geometry. The matter 
fields p, f and Q l are obtained after every step from 



A, 



(eukrj +eji k fi)f l J k 



(34) 



with f 1 the radial unit vector and J 1 the puncture's angu- 
lar momentum. Constructing initial data requires then 
solving only the Hamiltonian constraint fl9j 



1 



AVj + -^" 7 iiji y + 2 7n/> 5 p = 



(35) 



with Aij given by (|34|) and A the flat Laplacian. No- 
tice that the standard assumptions of conformal flatness 
and vanishing trace of the extrinsic curvature have been 
used. Once a solution to Eq. (|35|) is found, the initial 
data for the spatial metric and the extrinsic curvature 
are obtained from 7^ = ip 4 Vij an d K i3 = ijj~ 2 Aij, re- 
spectively. 

To solve Eq. (|35|) . we use the puncture ansatz ip = 
1 + to p /(2 r) + u. with m p the puncture's bare mass pa- 
rameter. We choose the source p as: 



p a e 



1 



2r 
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(36) 



and the bare angular momentum as J % = J z l . For all 
our simulations, we set J = 0.8. The last factor in ([36]) 
is needed for regularity at the location of the puncture. 
Thus, Eq. {35]) becomes: 
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J 2 r sin 2 6 



8 (r + ra p /2 + ur) 7 



-2tt 



r + m p /2 + ur 
r + m p /2 



p a e 



-(r-r a ) 2 /a 2 



0,(37) 



5 



with f l z 3 rjij — cos 6. 

The parameters r and a are chosen to favor the ac- 
cretion of most of the shell by the BH in evolution time- 
scales of < 40 M. The energy density amplitude parame- 
ter p Q is the knob that controls the mass associated with 
the shell, and thus regulates the total ADM mass M in 
the initial data. 

Our approach to breaking the extremality bound Q = 
J I (2 M H ) = 1 is to not only increase the angular momen- 
tum J of the BH, but most importantly to decrease its 
mass Mh- To accomplish this, we endow the shell sur- 
rounding the BH with a negative mass, which decreases 
the BH's mass as it gets accreted. 

It should be noted that in isolation any spinning punc- 
ture that models a rotating BH has a small (< 0.1% of 
the total energy) spurious amount of gravitational radi- 
ation. This spurious radiation carries away a burst of 
angular momentum t 20j ; this amount, however, is small 
and does not affect the conclusions from our numerical 
experiment. 



Run 


r /M 


a/M a 


PoMl 




m p /M 


Xo 


VO 











0.62 


0.80 


VI 


0.88 


0.88 


+7.03 x 10" 


3 


0.55 


0.62 


V2 


1.03 


1.03 


-3.77 x 10" 


4 


0.64 


0.85 


V3 


1.07 


1.07 


-1.05 x 10" 


3 


0.71 


0.90 


V4 


1.18 


1.18 


-1.72 x 10" 


3 


0.72 


1.12 


V5 


1.23 


1.23 


-2.00 x 10" 


3 


0.76 


1.20 



TABLE I: Initial parameters denning the matter shell 
(r ,a,p ). The last two columns display the puncture mass 
parameter m p and ADM spin parameter \o = Jo/M„. All 
parameters are given in units of the ADM mass M 



Table U shows the parameters used for our simulations: 
the shell parameters (r D , a, p a ) and the puncture mass 
m p . These quantities and further results are given in di- 
mcnsionless units in terms of the ADM mass M Q . The 
last column in Table Q] reports for each simulation the 
total spin parameter \ Q = J /M%, with J a the total 
ADM angular momentum. Since \o is associated with 
the ADM mass and angular momentum of the space- 
time, it is not subject to an extremality condition. In 
particular, Xo includes angular momentum from outside 
the horizon. Notice that case V0 consists of a single 
puncture in vacuum. This case is used as a control run. 
Case VI is a fiducial evolution with a positive energy den- 
sity shell. Cases V2-V5 contain shells with increasingly 
negative energy density. They are the central piece of 
our study. Table |H] gives the initial values for the irre- 
ducible Mh and Christodoulou M masses for each case. 
The table also shows the initial values of the BH's spin 
parameters \ = J/M 2 , ( = J/(2M H ), and its surface 
gravity n. 



Run 


M H /M 


M/M 


X 


c 


K M 


V0 


0.886 


0.994 


0.810 


0.511 


0.209 


VI 


0.794 


0.886 


0.796 


0.496 


0.237 


V2 


0.908 


1.021 


0.814 


0.515 


0.202 


V3 


0.958 


1.071 


0.800 


0.499 


0.196 


V4 


1.020 


1.156 


0.830 


0.533 


0.175 


V5 


1.056 


1.120 


0.835 


0.538 


0.168 



TABLE II: Initial irreducible mass Mh, Christodoulou mass 
M, dimensionless spin parameters \ = J/M 2 and £ = 
J/(2 M H ), and surface gravity k of the BH in units of the 
total ADM mass M a 



Although we are dealing with a single BH, the numer- 
ical simulations are challenging. As the BH approaches 
extremality, the horizon undergoes extreme pancake-like 
deformation. Because of this severe deformation, in or- 
der to capture the dynamics in the vicinity of the BH and 
to have a chance of locating its AH, we were forced to 
use meshes with 104 2 x 52 grid point shapes, and resolu- 
tions in the finest mesh of at least M/200, in addition to 
using 6th-order accurate finite differencing. For some of 
the cases, we carried out simulations with resolutions of 
M/167 in the finest mesh to investigate the dependence 
of our results with resolution. We did not find notice- 
able differences regarding the onset of super-extremality. 
Resolution effects were mostly visible in tracking the AH. 

V. SUPER-EXTREMAL BLACK HOLES 



0.15 




20 40 60 80 100 120 

t(M) 



FIG. 1: Fractional change of the Christodoulou mass (solid) 
and spin (dotted) for the cases of VI (top) and V4 (bottom) . 

We evolved the series of runs described in Table |U 



6 



As previously stated, the case VO is a control simula- 
tion to show that in the absence of matter, the BH mass 
and angular momentum indeed remain constant. For all 
the non-vacuum cases, the shell is absorbed early in the 
evolution. As a consequence, the mass Mh of the BH 
changes, increasing for the positive energy case VI and 
decreasing for the remaining negative energy cases V2- 
V5. 

The case VI involving a positive energy density shell 
demonstrates an important feature of our data: the shell 
of matter has negative angular momentum and positive 
mass (see Figure [IJ. The BH mass increases and the spin 
decreases as the shell is absorbed. (The time axis in all 
figures is given in units of the initial Christodoulou mass 
of the BH.) Figure [1] also shows the behavior of the BH 
mass and angular momentum for the V4 case with a neg- 
ative energy shell. The Equivalence Principle indicates 
that motion of test bodies is independent of their mass 
(even the sign of their mass). While our shells are not 
test bodies, we do expect the motion of the shells for the 
VI positive mass case and the V4 negative mass case to 
be similar. We thus anticipate and find, for example in 
V4, that as the negative mass shell falls into the BH, the 
BH mass decreases and the BH angular momentum in- 
creases (see Figure [T]), thus indicating that our negative 
energy shells have initially positive angular momentum. 




t(M) 



FIG. 2: Modes (2,0) and (2,2) of the Weyl scalar * 4 for all 
cases in the first 120 M of evolution time. 

Although the initial configuration is axisymmetric, 
there is in all cases gravitational radiation emitted during 
the evolution. The first emission occurs early on around 
t ~ 40 M. This is due to the spurious gravitational radi- 
ation in spinning punctures mentioned in the previous 
section. Details of this emission can be seen in Fig- 
ure [2] where we show the (2,0) and (2,2) modes of the 
Weyl scalar ^4. In addition, the accretion induces non- 
axisymmctric deformations that trigger a larger burst 
around t ~ 60 M, followed by quasi- normal ringing (see 




t(M) 



FIG. 3: Modes (2,0) and (2,2) of the Weyl scalar * 4 for all 
cases after 150 M of evolution time. 



top panel Fig. [5]) . Notice from the bottom panel of Fig. [2] 
that the (2,2) mode, which could potentially carry angu- 
lar momentum, is slightly stronger for the V4 and V5 
cases. At late times, t > 150 M, non- axisymmetric in- 
stabilities trigger a much stronger additional burst of ra- 
diation for those two cases, as seen in Figure |3] As a 
consequence, the BH looses angular momentum. 




100 150 200 250 300 350 400 450 

t(M) 

FIG. 4: Spin parameter £ as a function of evolution time. 

Figure |4] shows the core of our study, namely the evo- 
lution of the spin parameter £ during the spacetime evo- 
lution. As expected, C is constant in the V0 case. Also 
expected in case VI, with a positive energy density shell, 
is the decrease of the spin parameter ( of the BH due 
to the monotonic increase (decrease) of its mass (spin) 
observed in Fig. [TJ 

Cases V2-V5, with their negative energy density shells, 
are our main focus. As the shell is swallowed, the mass of 
the BH decreases monotonically since it accretes negative 
mass. Furthermore, given that our negative energy shells 
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have positive angular momentum, the spin of the BH will 
monotonically increase during the accretion. When taken 
together, the increase in spin and decrease in mass yield 
the observed growth of ( in Fig. @] Notice from Fig. @] 
that the larger the negative energy density of the shell, 
the faster the increase experienced by (. Also, only the 
cases V4 and V5 breach the extremality bound ( = 1. 
These cases are the ones in which Xo > 1 initially (see 
Table |T| . V4 becomes super-extremal at a time ~ 65 M 
and V5 at ~ 45 M. In both cases, C continues to grow un- 
til it reaches a maximum, at ~ 85 M for V4 and ~ 70 M 
for V5. The drop in £ after reaching the maximum in 
V4 and V5 is due to the emission of angular momentum. 
The BH is not able to sustain super-extremality and si- 
multaneously retain axisymmetry. 

Fig. 2] also shows that after V4 and V5 reach a maxi- 
mum, the £ lines terminate. This signals the time when 
we are no longer able to locate the AH of the BH. We 
have investigated whether the loss of the AH is real or 
due to numerical effects. Simulations at different resolu- 
tions tells us that most likely the loss of the horizon is a 
numerical artifact: the last horizon times depend on res- 
olution. The horizon undergoes a severe pancake defor- 
mation that the AH tracker is not able to handle because 
of the lack of resolution. We stress that the simulation 
does not crash; we only lose the horizon. This should not 
be viewed as "strange" since it is well known that stable 
puncture simulations do not resolve the puncture nor its 
immediate vicinity. If the horizon is too small, as with 
rapidly spinning BHs, the horizon is in danger of entering 
the under-resolved region near the puncture. 

Notice in Fig. [4] that for V4 around 150 M, when £ 
becomes again subcritical, we are again able to find the 
AH, albeit only briefly. Although the value of £ has sig- 
nificantly decreased, the shape of the horizon remains 
pancake-like, thus the challenge of locating the horizon 
remains. Since the AH shape is coordinate dependent, we 
are currently investigating modifications to the puncture 
gauges that alleviate the horizon deformations, giving us 
a better chance of locating the horizon. 

Finally, it is evident in Fig. |4] the increase that £ ex- 
periences in the sub-extremal cases V2 and V3. V2 gives 
hints of reaching an asymptotic value at later times. On 
the other hand, V3 has an approximate linear growth 
between 75 M and 300 M when it saturates at £ ~ 0.92. 



VI. NULL ENERGY CONDITION 

As mentioned in Sec. HH one of the assumptions that 
the extremality condition £ < 1 hinges on is the validity 
of the null energy condition. The null-energy condition 
states that for all null vectors k a 

N = T ab k a k b > . (38) 




FIG. 5: Snapshots of J\f computed from Eq. (|40)l . From top 
to bottom, case V3, V4 and V5 are shown. Left panels show 
a cut through the equatorial plane, and right panel along 
the rotation axis of the BH. White areas are regions with 
positive values of JV, and gray shaded areas those in which 
the null energy condition is violated. Snapshots are at times of 
65.28M, 62.72M, and 50.24M for V3, V4 and V5 respectively 
and cover a region 6M across. The grayscale is logarithmic 
in absolute value with a global minimum of -0.005. 

Substitution in Eq. (|3"5|) of the stress-energy tensor T a & 
with the form given by Eq. (|33p yields 

N = 2 -p-2j a l a + Ml F l a l b . (39) 

Above, we have used a null vector k a = n a + l a with 
l a n a = and l a l a = 1. If we choose l a to be radial and 
centered at the BH, i.e. l a = e", then 

N = 2 - P -2 ]r + Mj r F 7 (40) 

We have evaluated (|4"0)) throughout the computational 
domain. Not surprisingly, we found that for all the neg- 
ative energy density cases the null energy condition is 
violated in regions near the BH. Figure [5] shows for cases 
V3, V4 and V5 snapshots of both the equatorial plane 
(left panels) and the rotation axis (right panels) of the 
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FIG. 6: Average of the null energy condition Mh on the sur- 
face of the AH. 



BH. White areas are regions with positive values of TV, 
and gray shaded areas those in which the null energy 
condition is violated. The grayscale is logarithmic in ab- 
solute value, with a global minimum of -0.005. 

Notice in the left panels in Fig. [5] the evident inspi- 
ral structure in M associated with the accretion of the 
negative energy shell. Also interesting is that among the 
three cases depicted, only the one with the weakest neg- 
ative energy shell shows no regions where the null energy 
condition is not violated. 

Wc have also evaluated the null energy condition at the 
surface of the AH. That is, we evaluated Mh = T a bk a k b 
on the AH using the null vector k a = n a + s a , with s a 
the spatial unit normal to the AH. In Fig. E] we show 
how the average of Mh changes during the course of the 
evolutions. The null energy condition is clearly violated 
in cases V2-V5. A closer look at the data shows that this 
null energy condition is violated not only on average, but 
also everywhere on the AH surface. 

In summary, the violations of the null energy condition 
in the vicinity of the BH (Fig. [5]) and on the surface of the 
AH (Fig. demonstrates consistency with the violation 
of the extremality bound £ = 1. 



VII. CONSTRAINTS AND LATE BEHAVIOR 

To understand the late behavior and, in particular, to 
get clues about the instability of the super-extremal BHs, 
we calculated the L2 norms of p, j l and Q % within con- 
centric shells 1M < r < 2M, 2 M < r < AM, and 
4 M < r < 8 M as a function of evolution time. Fig- 
ures |5] and [5] show these L2 norms. First note that 
after t ~ 100 M the L2 norms of p, f and Q l for the 
V0, VI and V2 cases reach comparable levels. This sig- 
nals that the VI and V2 systems are settling down to a 
vacuum BH. Given that the BHs in VI and V2 accreted 
matter, their final mass and spins will be different from 
those of V0. 



Case V3 shows a different behavior. Initially, p, j z and 
Q l drop and stabilize as with VI and V2 as a consequence 
of matter having been accreted by the BH. The approxi- 
mately constant values are, however, larger than the cor- 
responding values in VI and V2. Beyond, t ~ 250 M, p, 
j l and Q % begin to grow as the BH approaches its maxi- 
mum £ value at t ~ 300 M (see Fig. @|. 



10" 1 
10" 2 




50 100 150 200 250 300 350 

t(M) 

FIG. 7: L2 Norm of p (r.h.s. Hamiltonian constraint) in 
concentric shells: 1 M < r < 2 M (solid), 2 M < r < 4M 
(dashed), and AM <r <8M (dotted). 

Cases V4 and V5, those that led to super-extremal 
BHs, show more complex behavior. Initially, the slight 
drop of p, j % and Q l in the outer two shells gives a hint of 
the undergoing accretion with the inner shell remaining 
approximately constant. Around the time the BH sur- 
passes extremality, t ~ 65 M for V4 and t ~ 45 M for 
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FIG. 8: L2 Norm of j* (r.h.s. momentum constraint) in 
concentric shells: 1M < r < 2M (solid), 2M < r < AM 
(dashed) , and 4 M < r < 8 M (dotted) . 
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FIG. 9: L2 Norm of the connection constraint Q a in concentric 
shells: 1M < r < 2M (solid), 2M <r < AM (dashed), and 
AM <r <8M (dotted). 



V5, p, f and Q l start growing appreciably. This smooth 
growth continues until t ~ 120 M for V4 and t ~ 100 M 
for V5. During this period, £ reaches a maximum value 
and proceeds to decrease. The decrease in £ is eventu- 
ally accompanied with a decrease in p, j % and Q l . At 
around t ~ 150 M for V4 and t - 130 M for V5, an 
ejection burst is triggered deep in the interior of the BH. 
Evidence of this ejected burst is the delay at which the 
burst emerges in each of the shells. That is, in Figures [71 
H] and |9] the burst shows in the solid line first, dash line 
next and dotted line last. 

The L2 norms in Figures [7J [5] and [S] do not provide 
a good sense of the details of the dynamics of p, j l 
and Q % . In Fig. [TU] we display the temporal evolution 



of spatial features. We show snapshots of p for the V4 
case in the equatorial plane (left) and axis plane (right) 
at times t ~ 50 M, 100 M, 150 M, and 200 M from top 
to bottom. The grayscale from top to bottom is such 
that [white:black] = [0 : -0.003], [0.08 : -0.08], [0.02 : 
-0.02], and [0.008 : -0.008] respectively. The panels are 
5 M across. The bottom left two panels clearly show a 
burst of negative energy. A similar burst is found in j % , 
which explains the smaller value of £ during a brief period 
at t ~ 150 M when the AH is again located. A detailed 
investigation of the late behavior of these super-extremal 
cases will be the subject of subsequent study. 
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FIG. 10: Snapshots of p of the V4 case in the equatorial-plane 
(left) and axis-plane (right) at times t ~ 50 M, 100 M, 150 M, 
and 200 M. The grayscale from top to bottom is such that 
[white:black] = [0 : -0.003], [0.08 : -0.08], [0.02 : -0.02], and 
[0.008 : —0.008] respectively. The regions, all centered about 
the coordinate origin, are 5 M across. 

VIII. CONCLUSIONS 

We have carried out a series of numerical experi- 
ments showing accreting BHs that violate the extremal 



spin condition £ < 1. The experiments consisted of a 
BH, modeled by a spinning puncture, surrounded by a 
spherically-symmetric shell. We considered shells with 
positive energy density and negative angular momentum, 
but the main focus was shells with negative energy den- 
sity and positive angular momentum. The idea behind 
this setup was that, as the BH accretes a negative energy, 
positive angular momentum shell, its mass will decrease 
and its spin will increase, leading to possible violations 
of the C = 1 extremality bound. We were successful in 
violating the extremality bound, at least temporarily. In 
agreement with the findings by Booth and Fairhurst , 
the violations were accompanied with violations of the 
null energy condition. The BHs that violated the £ = 1 
bound were not able to sustain super-extremality and si- 
multaneously retain axisymmctry. 

The study involved several challenges. The most sig- 
nificant challenge was locating the AH as the BH became 
super-extremal. In subsequent work, we will investigate 
gauge conditions that could alleviate difficulties caused 
by the extreme deformations of the BH horizon. If a new 
gauge condition is found, we will be in a better position to 
investigate if violations of extremality are accompanied 
with the disappearance of the AH. We will also focus on 
the late-time behavior to get a better understanding of 
the causes behind the late emission responsible for the 
drop of the spin parameter £. 

Our study introduced a new approach to construct 
non-vacuum dynamical spacctimes: the matter-without- 
matter evolution framework. Under this approach, only 
the evolution equations are used, with the matter source 
terms set to zero. Specifically, we demonstrated that the 
3+1, BSSN vacuum evolution equations were capable of 
providing the dynamics of matter fields that are invisible 
to the equations. The price paid in MWM evolutions is 
restrictions on the "equations of state" satisfied by the 
matter content. It is not clear whether MWM evolu- 
tions could be applied to other 3+1 formulations of the 
Einstein equations and produce stable and convergent 
simulations. 
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